Josephson effect in graphene SBS junctions 
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We study Josephson effect in graphene superconductor- barrier- superconductor junctions with 
short and wide barriers of thickness d and width L, which can be created by applying a gate voltage 
Vb across the barrier region. We show that Josephson current in such graphene junctions, in complete 
contrast to their conventional counterparts, is an oscillatory function of both the barrier width d and 
the applied gate voltage Vo- We also demonstrate that in the thin barrier limit, where Vb — ► oo and 
d — ► keeping Vod finite, such an oscillatory behavior can be understood in terms of transmission 
resonance of Dirac-Bogoliubov-de Gennes quasiparticles in superconducting graphene. We discuss 
experimental relevance of our work. 
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I. INTRODUCTION 

Graphene, a two-dimensional single layer of graphite, 
has been recently fabricated by Novoselov et. alX. In 
graphene, the energy bands touch the Fermi energy at 
six discrete points at the edges of the hexagonal Bril- 
louin zone. Two of these six Fermi points, referred to 
as if and K' points, are inequivalent and the quasiparti- 
cle excitations about them obey linear Dirac-like energy 
dispersion 2 . The presence of such Dirac-like quasipar- 
ticles leads to a number of unusual electronic proper- 
ties in graphene including relativistic quantum hall ef- 
fect with unusual structure of Hall plateaus 3 -, which has 
been verified in experiments 4 . Further, as suggested in 
Ref. [H, Dirac quasiparticles in graphene leads to real- 
ization of interesting physical phenomenon such as Klein 
parados* 4 ^, Lorenz-boost type phenomenon 7 -, and uncon- 
ventional Kondo effect 8 -^. 

Another interesting consequences of the existence 
Dirac-like quasiparticles can be understood by studying 
superconductivity in graphene. It has been suggested 
that superconductivity can be induced in a graphene 
layer in the presence of a superconducting electrode 
near it via proximity effec t 10 ! 11 ' 12 or by intercalating it 
with metallic atoms 13 . Consequently, studies on tunnel- 
ing conductance on both normal metal-superconductors 
(NS) and normal metal-barrier-superconductor (NBS) 
junctions in graphene have been undert ake n 1 1 14 1 1 5 . It 
has been shown in Refs. [U and [TH that the tunneling 
conductance of such NBS junctions are oscillatory func- 
tions of the effective barrier strength and that this oscil- 
latory phenomenon can be understood in terms of trans- 
mission resonance phenomenon of Dirac-Bogoliubov-de 
Gennes (DBdG) quasiparticles of graphene. Josephson 
effect has also been studied in a superconductor-normal 
metal-superconducting (SNS) junction in graphen e 12 ' 16 . 
It has been shown in Ref. [l2], that the behavior of such 
junctions is similar to conventional SNS junctions with 
disordered normal region. Such Josephson junctions with 
thin barrier regions have also been experimentally real- 
ized recentlyi 7 -. However, Josephson effect in graphene 
superconductor-barrier-superconductor (SBS) junctions 



has not been studied so far. 

In this work, we study Josephson effect in graphene 
for tunnel SBS junctions. In this study, we shall con- 
centrate on SBS junctions with barrier thickness d -C £ 
where £ is the superconducting coherence length, and 
width L which has an applied gate voltage Vb across the 
barrier region. Our central result is that in complete 
contrast to the conventional Josephson tunnel junctions 
studied so fa r 18 ' 19 , the Josephson current in graphene 
SBS tunnel junctions is an oscillatory function of both 
the barrier thickness d and the applied gate voltage Vq. 
We provide an analytical expression for the Josephson 
current of such a junction. We also compute the crit- 
ical current of graphene SBS junctions. We find that 
this critical current is also an oscillatory function of Vb 
and d and study the amplitude and periodicity of its os- 
cillation. We also show that in the thin barrier limit, 
where the barrier region can be characterized by an ef- 
fective dimensionless barrier strength \ — Vod/hvF (u.f 
being the Fermi velocity of electrons in graphene), the 
Josephson current becomes an oscillatory function of \ 
with period it— . We find that in this limit, the oscil- 
latory behavior of Josephson current can be understood 
as a consequence of transmission resonance phenomenon 
of Dirac-Bogoliubov-de Gennes (DBdG) quasiparticles in 
graphene. We demonstrate that the Josephson current 
reaches the Kulik-Omelyanchuk limit 20 for \ — nir (n be- 
ing an integer), but, unlike conventional junctions, never 
reaches the Ambegaokar-Baratoff limit^l for large X- We 
also discuss simple experiments to test our theory. 

The organization of the rest of the paper is as follows. 
In Sec. HH we obtain an analytical expression for Joseph- 
son current for a general SBS junction of thickness d -C £ 
and applied voltage Vb and demonstrate that the Joseph- 
son current is an oscillatory function of both d and Vq. 
This is followed by Sec. lIIIi where we discuss the limiting 
case of thin barrier and demonstrate that the oscillatory 
behavior of the Josephson current can be understood in 
terms of transmission resonance of DBdG quasiparticles 
in graphene. Finally we discuss experimental relevance 
of our results in Sec. HVl 
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II. JOSEPHSON CURRENT FOR TUNNEL SBS 
JUNCTIONS 

We consider a SBS junction in a graphene sheet of 
width L lying in the xy plane with the superconduct- 
ing regions extending x = -oo to j; = — d and from 
x = x — oo to for all < y < L, as shown in Fig. Q] 
The superconducting regions x > and x < —d shall be 
assumed to be kept close to superconducting electrodes 
so that superconductivity is induced in these regions via 
proximity effec t 10 ! 11 . Alternatively, one can also possibly 
use intercalated graphene which may have s-wave super- 
conducting phases^. In the rest of this work, we shall 
assume that these regions are superconducting without 
worrying about the details of the mechanism used to in- 
duce superconductivity. The region B, modeled by a 
barrier potential Vo, extends from x = — d to x = 0. 
Such a local barrier can be implemented by either us- 
ing the electric field effect or local chemical doping^. 
In the rest of this work, we shall assume that the bar- 
rier region has sharp edges on both sides which requires 
d « A = 2ir /hp, where kp and A are the Fermi wave- 
vector and Fermi wavelength for graphene. Such barriers 
can be realistically created in experiments^. The width 
L of the sample shall be assumed to be large compared to 
all other length scales in the problem. The SBS junction 
can then be described by the DBdG equations^ - 

f-H a -E F + U(r) A(r) \ 

^ A*(r) E F - U{r) - H a ) Va ~ Va ' 

(1) 

Here , ^„ = „ , a , ^ s , a ) are the 4 component 
wavefunctions for the electron and hole spinors, the in- 
dex a denote K or K' for electron/holes near K and K' 
points, a takes values K'iK) for a = K(K'), Ep denote 
the Fermi energy, A and B denote the two inequivalent 
sites in the hexagonal lattice of graphene, and the Hamil- 
tonian TL a is given by 

H a = -ihv F (a x d x + sgn(a)a y d v ) , (2) 

where sgn(a) takes values ± for a = K(K'). 

The pair-potentials A(r) in Eq.[T]connects the electron 
and the hole spinors of opposite Dirac points. We have 
modeled the pair-potential as 

A(r) = A [exp(# 2 )0(» + exp(ifa)0(x + d)] , (3) 

where Ao is the amplitude and <pU2) are the phases of 
the induced superconducting order parameters in regions 
I (II) as shown in Fig. [2 and 9 is the Heaviside step func- 
tion. Notice that the mean-field conditions for supercon- 
ductivity is satisfied as long as Ao -C Ep or equivalently 
kpt; 3> 1, where £ = Hvp/ttAq is the superconducting 
coherence length^. The potential U (r) gives the relative 
shift of Fermi energies in the barrier and superconducting 
regions and is modeled as 

U(r) = V Q 6(-x)6(x + d). (4) 
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FIG. 1: A schematic graphene SBS junction with the barrier 
B sandwiched between two superconductors I and II with pair 
potentials Aoe^ 1 and A e i</>2 . The barrier region is created 
by an external gate voltage Vo. 

Solving Eq. [I] we obtain the wavefunctions in the su- 
perconducting and the barriers regions. In region I, for 
the DBdG quasiparticles moving along ±x direction with 
a transverse momentum k y = q = 2itn/L (for integer n) 
and energy e, the wavefunctions are given byi°- 

^ = (uf, u± uf, «±)e«^-*+w>+«x], (5 ) 

where 

-| = ±exp(±i 7 ), -|=exp[-i(<£iT/3)], 
u ± 

-4 = ±exp[±z( T fth+/3 + 7)], (6) 
u x 

and X)i=i4l u i| 2 — is the normalization condition 
for the wavefunction for d -C k _1 , where k _1 = 
(HvF) 2 k s / [E F Aq sin(/3)] is the localization length. Here 

k s = *J (Ep/hvF) 2 — q 2 , 7, the angle of incidence for the 
quasiparticles, is given by sin(7) = hvpq/Ep, and f3 is 
given by 

P = coe-^e/Ao) if|e|< A , 

= -icoeh-^e/Ao) if |e| > A , (7) 

Note that for |e| > Ao, K becomes imaginary and the 
quasiparticles can propagate in the bulk of the supercon- 
ductor. The wavefunctions in region II (x > ) can also 
be obtained in a similar manner 

i>fi = (vf,vi,vi,vf)eM ±k «*+^-™\ (8) 

where Sj=i 4 l^l 2 = %k and the coefficients Vi are given 
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by 



± exp(±ry), -|- = exp[-i(0 2 ± /?)], 



= ±exp[±i(=F0i -/9 + 7)]> 



(9) 



The wavefunctions for electrons and holes moving 
along ±x in the barrier region is given by 



e± 



1, ±e ±ie , 0, 0) exp [i (±fe b a; + qy)] /V2d, 



^ = (o,0,l,±e^ e ')exp[i(±k' b x + qy)}/V2d. 



(10) 

Here the angle of incidence of the electron(hole) 0(9') is 
given by 



sin [9(0')] 



h(K) = 



hvpq 



e+(-)(E F -V Q ) 



'e+(-)(E F -V ) 



-q 2 (11) 



To compute the Josephson current in the SBS junc- 
tion, we now find the energy dispersion of the subgap 
Andreev bound states which are localized with localiza- 
tion length at the barrie r 22 ' 23 . The energy dispersion 
e„ (corresponding to the subgap state characterized by 
the quantum number n) of these states depends on the 
phase difference <f> = 4> 2 — 4>i between the superconduc- 
tors. It is well known that the Josephson current / across 
the junction at a temperature T is given b y 12 ' 22 



!{<P;x,T ) 



4e 
1 



E E 



de n 



f(tn), (12) 



where f(x) = l/(e x '( kBT °' + 1) is the Fermi distribution 
function and ks is the Boltzman constant^. 

To obtain these subgap Andreev bound states, we now 
impose the boundary conditions at the barrier. The 
wavefunctions in the superconducting and barrier regions 
can be constructed using Eqs. [5l [8l and [TOl as 



aiipf+biipj vp// = a 2 tpjj + b 2 ijj 



II- 



*j 

* B = w4 + + #r + + #£r> (13) 

where 01(02) and 61(62) are the amplitudes of right and 
left moving DBdG quasiparticles in region I(II) and p(q) 
and r(s) are the amplitudes of right(left) moving electron 
and holes respectively in the barrier. These wavefunc- 
tions must satisfy the boundary conditions: 



=-di 



B x=0 



II x=0- 



(14) 



Notice that these boundary conditions, in contrast their 
counterparts in standard SBS interfaces^ 3 -, do not impose 
any constraint on derivative of the wavefunctions. Thus 



the standard delta function potential approximation for 
short barrier a 22 ' 23 can not be taken the outset, but has 
to be taken at the end of the calculations. 

Substituting Eqs. [5J E EE and [JJ in Eq. QH we find 
the equations for boundary conditions at x = — d to be 



die 



-ik s d—Kd , i ik s d—Kd 



6ie lfes 



p e - ik bd _|_ q e ik bd 



i(y — k s d) — Kd ^^ — i{ r ) — k s d) — Kd 

pe t(e-k b d) _ -i(e-k b d) 



aie -i(4>l-l3+k s d)-Kd _|_ i(0i+|3-fc s d)-red 



l ^ e i(l-<t>l+P-k 3 d)-nd _ j } ^ e -i(~f+4> 1 +f3-k B d)-Kd 



re -i(6'+k b d) _ se i{6'+k' b d)_ / 15 ^ 



Similar equations at x = reads 

02 + 62 = p + q 

a2e n _ b 2 e~^ = pe w - qe- %0 

a 2e -*(*»+«+6 2e -*C*»-« = r + s 

a2e *(7-02-/3) _ &2e -i(7+^2-« = re -i0' _ se i0' ( 16 ) 

Eqs . fl5l and [TBI therefore yields eight linear homogeneous 
equations for the coefficients 0^=1,2, 6i=i,2 5 9i r i and 
s, so that the condition for non-zero solutions of these 
coefficients can be obtained as 



A' sin(2/3) + B' cos(2/3) + C = 
where A', B', andC are given by 



(17) 



A' 



B' = 



C = 



cos 

+ COS 



(k' b d) cos(7) cos(0') sin(fcf,d) (sin(7) sin(0) — 1) 
Ms(kbd) cos(7) cos(0) sin(k' b d) 

+ — cos(ki)d) cos(0) sin(27) sin(0') sm(k' b d) 

sm(k' b d) sin(fcfed) [ — 1 4- sin(0) sin(7) 

- sin(0') sin(7) + sin(0) sin(0') sin 2 (7)] 

— cos(k b d) cos(k' b d) cos 2 (7) cos(0) cos(0') 
cos 2 (7) cos(0) cos(0') cos(</)) — sm(kbd) sin(k' b d) 
x [sin(0) shi(0') - sin 2 (7) 

-sin( 7 ) (sin(0) - sin(0'))] 



(18) 



Note that in general the coefficients A', B', and C de- 
pends on e through kb, k' b , and 0' which makes it 
impossible to find an analytical solution for Eq. 1171 
However, for subgap states in graphene SBS junctions, 
e < Ao <C Ep. Further, for short tunnel barrier we 
have \Vq — Ep\ > Ep- In this regime, as can be seen 
from Eqs. [TTJ A', B' , and C' become independen t of 
e since kb — k' b ~ k\ = \/[(Ep — Vo)/hvp} 2 — q 2 and 
~ -0' ~ 0i = sin -1 [hv F q/(E F - V )} so that the e 
dependence of k b , k' b , and 0' can be neglected. In this 
regime one finds that A' ,B',C — > A,B,C where 



4 




FIG. 2: Plot of Josephson current 7 as a function of phase dif- 
ference 4> and the applied gate voltage Vo for k B T = 0.01A FIG. 3: Plot of I c /I vs the applied gate voltage Vo and the 
and d = 0.5A showing oscillatory behavior of I/I as a func- junction thickness d for T = 0.01 A . 
tion of the applied gate voltage. 



A 
B 



sin 2 (fcid) [1 — sin(7) sin(#i)] z 
cos 2 (kid) cos 2 (7) cos 2 (0i) 



C = sin 2 (fcid)[sin(7)-sin(6' 1 )] 1 
+ cos 2 (7) cos 2 (Q\ ) cos(</>) 



(19) 



The dispersion of the Andreev subgap states can now be 
obtained from Eqs. [T7J and [JJ One finds that there are 
two Andreev subgap states with energies e± = ±e where 



e = A av / l/2-C/2B 



(20) 



Using Eq. 1121 one can now obtain the expression for the 
Josephson current 

I(<p,V ,d,T Q ) = I g(<p,V ,d,T a ), 

cos 3 (7) cos 2 (0i ) sin(0) 



g((p,V ,d,T ) 



<~h 



Be/A 



x tanh(e/2fc B T ) 



(21) 



where Iq = eAQEpL/2h 2 irvF and we have replaced 
J2 q — ¥ EfL/(2ttHvf) J^, 2 dj cos(7) as appropriate for 
wide junctions^ 2 -. 

Eqs. [20l and [21] represent the central result of this 
work. From these equations, we find that both the dis- 
persion of the Andreev subgap states and the Joseph- 
son current in graphcnc SBS junctions, in complete con- 
trast to their conventional counterpart s 18 ! 19 ' 22 , is oscil- 
latory function of the applied gate voltage V and the 
barrier thickness d. This statement can be most easily 



checked by plotting the Josephson current I as a func- 
tion of the phase difference <p and the applied gate volt- 
age Vo for a representative barrier thickness d = 0.5A 
and temperature fcsT = 0.01A , as done in Fig. [2j In 
Fig. [31 we plot the critical current of these junctions 
I c (V ,d,T ) = Meac{I((f>,Vo,d,To)] as a function of the 
applied gate voltage Vq and barrier thickness d for low 
temperature UbTq = O.OlAo- We find that the critical 
current of these graphene SBS junctions is an oscillatory 
function of both Vo and d. This behavior is to be con- 
trasted with those of conventional junctions where the 
critical current is a monotonically decreasing function 
of both applied bias voltage Vo and junction thickness 

^18.19.22 

Next, we analyze the temperature dependence of the 
amplitude of oscillations of I c . To find the amplitude of 
oscillation, we have computed I c as a function of Vq (for 
a representative value of d = 0.3A), noted the maximum 
(7™ ax ) and minimum (/™ ln ) values of I c , and calculated 
the amplitude 7™ ax — 7™ m . The procedure is repeated for 
several temperatures T and the result is plotted in Fig. 0] 
which shows that the amplitude of oscillations decreases 
monotonically as a function of temperature. 

Finally, we discuss the period of oscillation of the crit- 
ical current. To obtain the period, we obtain the critical 
current I c as a function of barrier width d for the fixed 
applied gate voltage Vq and note down d pcr iod- We then 
compute Xperiod = V Q d pcl . iod /hv F and plot x pC riod as a 
function of Vo for fc^To = O.OlAo as shown in Fig. [5] 
We find that Xperiod decreases with Vq and approaches 
an universal value 7r for large Vo > 20Ep. This property, 
as we shall see in the next section, can be understood 
by analysis of graphene SBS junctions in the thin bar- 
rier limit (Vb — > 00 and d — > such that x — Vad/hvF 
remains finite^) and is a direct consequence of transmis- 
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0.4 0.6 

k B T /A 

FIG. 4: Plot of the temperature dependence of the amplitude 
of oscillations of I c (given by [J c max (d) - I c min (d)]// ) for d = 
0.3A. The amplitude is measured by noting the maximum 
and minimum values of the critical current by varying Vb for 
a fixed d. 



sion resonance phenomenon of DBdG quasiparticlcs in 
superconducting graphene. 



III. THIN BARRIER LIMIT 

In the limit of thin barrier, where Vq — ► oo and d — ► 
such that x — Vod/hvF remains finite, 9\ — > and k\d — > 
X- From Eqs. [T9l and [20l we find that in this limit, the 
dispersion of the Andreev bound states becomes 



4 b (<^;x) = ±Ao^l-T( 7 ,x)sin 2 (>/2), (22) 
cos 2 (7) 



r(7.x) 



1 — cos 2 (x) sin 2 (7) 



(23) 



where the superscript 'tb' denote thin barrier limit. The 
Josephson current I can be obtained substituting Eq. [221 
in Eq. [T2] In the limit of wide junctions, one gets 



9 th (^X,To) 



W b (0,X,T o ), 



tt/2 



d", 



T( 7 , x) cos( 7 ) sin(0) 



*/a L^/l-r( 7 , X ) sin 2 (0/2) 



tanh(e + /2fc B T ) 



(24) 



Eqs.[23l and [24l represent the key results of this section. 
From these equations, we find that the Josephson 
current in graphene SBS junctions is a 1 periodic 
oscillatory function of the effective barrier strength x 
in the thin barrier limit. Further we observe that the 
transmission probability of the DBdG quasiparticlcs in 
a thin SBS junction is given by T( 7 , x) which is also 
the transmission probability of a Dirac quasiparticle 
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FIG. 5: Plot of Xperiod of the critical current I c as a function 
of Vb- Note that Xperiod approaches n as we approach the thin 
barrier limit. 



through a square potential barrier as noted in Ref. [5]. 
Note that the transmission becomes unity for normal 
incidence (7 = 0) and when x = n7r - The former 
condition is a manifestation of the Klein paradox 
for DBdG quasiparticles^. However, this property is 
not reflected in the Josephson current which receives 
contribution from quasiparticles approaching the junc- 
tion at all angles of incidence. The latter condition 
(x = rm) represents transmission resonance condition 
of the DBdG quasiparticlcs. Thus the barrier becomes 
completely transparent to the approaching quasipar- 
ticles when x = nn an d in this limit the Josephson 
current reduces to its value for conventional tunnel junc- 
tions in the Kulik-Omelyanchuk limit: I th (<fi, nir, Tq) — 
4/ sin(0/2)Sgn(cos(0/2)) tanh (A |cos(0/2)| /2k B T )22.. 
This yields the critical Josephson current I* b (x = n7T ) — 
4Jo for UbTq <C Ao- Note, however, that in contrast to 
conventional junctions T( 7 , x) can n °t be made arbi- 
trarily small for all 7 by increasing \. Hence 7* b never 
reaches the Ambegaokar-Baratoff limit of conventional 
tunnel junctions 2 ^. Instead, /' b (x) becomes a ir periodic 
oscillatory function of x- The amplitude of these 
oscillations decreases monotonically with temperature 
as discussed in Sec. HU 

Finally, we compute the product I^Rn which is 
routinely used to characterize Josephson junction a 18 ' 19 , 
where Rn is the normal state resistance of the junction. 
For graphene SBS junctions Rn corresponds to the resis- 
tance of a Dirac quasiparticle as it moves across a normal 
metal-barrier-normal metal junction. For short and wide 
junctions discussed here, it is given by Rn = i?o/si(x) 
where R Q = -k 2 vf^ 2 / (e 2 E F L) and si(x) is given by^ii 2 - 



si(x) 



Tl/2 



tt/2 



d7 r (7,x)-COs(7). 



(25) 



Note that si(x) and hence Rn is an oscillatory func- 
tion of x with minimum 0.5i?o at x — n7T an d maximum 
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1.4-1 , , , , , , , , , , 

2 4 6 8 10 

X 

FIG. 6: Plot of ll h R N as a function of X - I^Rn is 
an oscillatory bounded function of x an d never reaches its 
value (7rAo/2e) for conventional junctions in the Ambegaokar- 
Baratoff limit. 



0.75i? at x = (n + V 2 ) 71 "- Tb -e product I^Rn, for thin 
SBS junctions is given by 

I?R N = (nA /2e)g^ ax ( X ,T)/ Sl ( X ), (26) 

where g^axix) is the maximum value of g th (4>, x). Note 
that /' b i?7v is independent of Ep and hence survives in 
the limit E F -> 12 . For k B T < A , ^(mr) = 4 
and si(mr) = 2, so that I^RN\x=nn = TtAo/e which 
coincides with Kulik-Omelyanchuk limit for conventional 
tunnel junction o 20 ' 23 . However, in contrast to the con- 
ventional junction, /* b i?jv for graphene SBS junctions do 
not monotonically decrease to the Ambegaokar-Baratoff 
limi t 21 i 23 of 7rAo/2e ~ 1.57Ao/e as \ is increased, but 



demonstrates 7r periodic oscillatory behavior and remains 
bounded between the values 7rA /e at X — n7T an d 
2.27A /e at x = {n + l/2)n, as shown in Fig. [6] 



IV. EXPERIMENTS 

As a test of our predictions, we suggest measuring 
DC Josephson current in these junctions as a function 
of the applied voltage Vq. Such experiments for con- 
ventional Josephson junctions are well-known 25 . Further 
SNS junctions in graphene has also been recently been ex- 
perimentally created^. For experiments with graphene 
junctions which we suggest, the local barrier can be fab- 
ricated by applying an additional gate voltage in the 
normal region of the junctions studied in Ref. ITtI In 
graphene, typical Fermi energy can reach Ep < 80meV 
with Fermi- wavelength A = 2n/kp > lOOnm 5 -. Effective 
barrier strengths of 500 — lOOOmeV and barrier widths of 
d ~ 20— 50nm can be achieved in realistic experiments^. 
These junctions therefore meet our theoretical criteria: 
d <C A and | Vb — Ep\ > Ep. To observe the oscillatory 
behavior of the Josephson current, it would be necessary 
to change Vb in small steps SVo. For barriers with fixed 
d/X = 0.3 and Vo/Ep = 10, this would require chang- 
ing Vq in steps of approximately 30meV which is exper- 
imentally feasible. The Joule heating in such junctions, 
proportional to I^Rn, should also show measurable os- 
cillatory behavior as a function of Vq. 

In conclusion, we have shown that the Josephson cur- 
rent in graphene SBS junction shows novel oscillatory 
behavior as a function of the applied bias voltage Vq and 
the barrier thickness d. In the thin barrier limit, such a 
behavior is the manifestation of transmission resonance of 
DBdG quasiparticles in superconducting graphene. We 
have suggested experiments to test our predictions. 

KS thanks V. M. Yakovenko for discussions. 
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